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Abstract 

We provide new evidence that the holographic dual to a strongly coupled charged Fermi liquid has 
a non-zero fermion density in the bulk. We show that the pole-strength of the stable quasiparticle 
characterizing the Fermi surface is encoded in the spatially averaged AdS probability density of 
a single normahzable fermion wavefunction in AdS. RecaUing Migdal's theorem which relates the 
pole strength to the Fermi-Dirac characteristic discontinuity in the number density at up^ we 
conclude that the AdS dual of a Fermi liquid is described by occupied on-shell fermionic modes in 
AdS. Encoding the occupied levels in the total probability density of the fermion field directly, we 
show that an AdS Reissner-Nordstrom black hole in a theory with charged fermions has a critical 
temperature, at which the system undergoes a first-order transition to a black hole with a non- 
vanishing profile for the bulk fermion field. Thermodynamics and spectral analysis confirm that the 
solution with non-zero AdS fermion-profile is the preferred ground state at low temperatures. 
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1 Introduction 



Fermionic quantum criticality is thought to be an essential ingredient in the full theory of high Tc 
superconductivity [I1I2]- The cleanest experimental examples of quantum criticality occur in heavy- 
fermion systems rather than high cuprates, but the experimental measurements in heavy fermions 
raise equally confounding theoretical puzzles [3]. Most tellingly, the resistivity scales linearly with 
the temperature from the onset of superconductivity up to the crystal melting temperature ^ and 
this linear scaling is in conflict with single correlation length scaling at criticality [5]. The failure 
of standard perturbative theoretical methods to describe such behavior is thought to indicate that 
the underlying quantum critical system is strongly coupled [SIE!- 

The combination of strong coupling and scale-invariant critical dynamics makes these systems 
an ideal arena for the application of the AdS/CFT correspondence: the well-established relation 
between strongly coupled conformal field theories (CFT) and gravitational theories in anti-de Sitter 
(AdS) spacetimes. An AdS/CFT computation of single-fermion spectral functions — which are 
directly experimentally accessible via Angle- Resolved Photoemission Spectroscopy lElElIin] — bears 
out this promise of addressing fermionic quantum criticality [HI [121 [131 El] (see also [ISIIIS])- The 
AdS/CFT single fermion spectral function exhibits distinct sharp quasiparticle peaks, associated 
with the formation of a Fermi surface, emerging from a scale-free state. The fermion liquid which this 
Fermi surface captures is generically singular: it has either a non-linear dispersion or non-quadratic 
pole strength pTl [TH] . The precise details depend on the parameters of the AdS model. 

From the AdS gravity perspective, peaks with linear dispersion correspond to the existence of 
a stable charged fermionic quasinormal mode in the spectrum of a charged AdS black hole. The 
existence of a stable charged bosonic quasinormal mode is known to signal the onset of an instability 
towards a new ground state with a pervading Bose condensate extending from the charged black 
hole horizon to the boundary of AdS. The dual CFT description of this charged condensate is 
spontaneous symmetry breaking as in a superfiuid and a conventional superconductor pT f [TSIIT^ [20] . 
For fermionic systems empirically the equivalent robust T = groundstate is the Landau Fermi 
Liquid — the quantum groundstate of a system with a finite number of fermions. The existence 
of a stable fermionic quasinormal mode suggests that an AdS dual of a finite fermion density state 
exists. 

We construct here the AdS/CFT rules for CFTs with a finite fermion density. The essential 
ingredient will be the Migdal's theorem, which relates the characteristic jump in fermion occupation 
number at the energy up of the highest occupied state to the pole strength of the quasiparticle. 
The latter we know from the spectral function analysis and its AdS formulation is therefore known. 
Using this, we show that the fermion number discontinuity is encoded in the probability density of 
the normalizable wavefunction of the dual AdS fermion field. 

These AdS/CFT rules prove that the AdS dual of a Fermi liquid is given by a system with 
occupied fermionic states in the bulk. The Fermi liquid is clearly not a scale invariant state, but 
any such states will have energy, momentum/pressure and charge and will change the interior 
geometry from AdS to something else. Which particular (set of) state(s) is the right one, it does 
not yet tell us, as this conclusion relies only on the asymptotic behavior of fermion fields near the 
AdS boundary. Here we shall take the simplest such state: a single fermion]^ Constructing the 
associated backreacted asymptotically AdS solution, we find that it is already good enough to solve 
several problems of principle: 

^ These solutions are therefore the AdS extensions of [HI [221 [231 [13 • 
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• Fermionic quantum critical systems should undergo a phase transition to a Fermi liquid at 
low temperatures, except when one is directly above the QCP. Likewise we find that the dual 
of the quantum critical state, a charged AdS black hole in the presence of charged fermionic 
modes, has a critical temperature below which fermionic Dirac "hair" forms. The derivative 
of the free energy has the characteristic discontinuity of a first order transition. This has to 
be the case: A fermionic quasinormal mode can never cause a linear instability indicative of 
a continuous phase transition. In the language of spectral functions, the pole of the retarded 
Green's function can never cross to the upper-half plane [I3]j^ The absence of a perturbative 
instability between this conjectured Dirac "black hole hair" solution and the "bald" charged 
AdS black hole can be explained if the transition is a first order gas-liquid transition. The 
existence of first order transition follows from a thermodynamic analysis of the free energy 
rather than a spectral analysis of small fiuctuations. 

• This solution with finite fermion number is the preferred ground state at low temperatures. 
The bare charged AdS black hole in a theory with charged fermions is therefore a false vacuum. 
Confusing a false vacuum with the true ground state can lead to anomalous results. Indeed 
the finite temperature behavior of fermion spectral functions in AdS Reissner-Nordstrom, 
exhibited in the combination of the results of [TTl [T3] and [12] , shows strange behavior. The 
former [HI [TJ] found sharp quasiparticle peaks at a frequency up = in natural AdS units, 
whereas the latter [12j found sharp quasiparticle peaks at finite Fermi energy cop 7^ 0. As 
we will show, both peaks in fact describe the same physics: the up 7^ peak is a finite 
temperature manifestation of (one of the) u = peaks in [13]. Its shift in location at finite 
temperature is explained by the existence of the nearby true finite fermion density ground 
state, separated by a potential barrier from the AdS Reissner-Nordstrom solution. 

• The charged AdS-black hole solution corresponds to a CFT system in a state with large ground 
state entropy. This is the area of the extremal black- hole horizon at T = 0. Systems with large 
ground-state entropy are notoriously unstable to collapse to a low-entropy state, usually by 
spontaneous symmetry breaking. In a fermionic system it should be the collapse to the Fermi 
liquid. The final state will generically be a geometry that asymptotes to Lifschitz type, i.e. 
the background breaks Lorentz-invariance and has a double-pole horizon with vanishing area, 
as expounded in [25] . Although the solution we construct here only considers the backreaction 
on the electrostatic potential, we show that the gravitational energy density diverges at the 
horizon in a similar way as other systems that are known to gravitationally backreact to a 
Lifshitz solution. The fully backreacted geometry includes important separate physical aspects 
— it is relevant to the stability of the Fermi liquid — and will be considered in a companion 
article. 

The Dirac hair solution thus captures the physics one expects of the dual of a Fermi liquid. We 
have based its construction on a derived set of AdS/CFT rules to describe systems at finite fermion 
density. Qualitatively the result is as expected: that one also needs occupied fermionic states in 
the bulk. Knowing this, another simple candidate is the dual of the backreacted AdS-Fermi-gas 

/electron star [26j which appeared during the course of this workj^ The difference between the 
two approaches are the assumptions used to reduce the interacting Fermi system to a tractable 
solution. As explained in the recent article [30], the Fermi-gas and the single Dirac field are the 

^Rcf. [39j argues that the instability can be second order. 

•^See also [27l[28]. An alternative approach to back-reacting fermions is |29) . 
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two "local" approximations to the generic non-local multiple fermion system in the bulk, in very 
different regimes of applicability. The electron-star/Fermi-gas is considered in the Thomas- Fermi 
limit where the microscopic charge of the constituent fermions is sent to zero keeping the overall 
charge fixed, whereas the single Dirac field clearly is the 'limit' where the microscopic charge equals 
the total charge in the system. This is directly evident in the spectral functions of both systems. 
The results presented here show that each pole in the CFT spectral function corresponds to a unique 
occupied Fermi state in the bulk; the electron star spectra show a parametrically large number of 
poles [311 |32l [30] , whereas the Dirac hair state has a single quasiparticle pole by construction. The 
AdS-Dirac-hair black hole derived here therefore has the benefit of a direct connection with a unique 
Fermi liquid state in the CFT. This is in fact the starting point of our derivation. 

In the larger context, the existence of both the Dirac hair and backreacted Fermi gas solution 
is not a surprise. It is a manifestation of universal physics in the presence of charged AdS black 
holes. The results here, and those of [HI [131 ESI [26] , together with the by now extensive literature 
on holographic superconductors, i.e. Bose condensates, show that at sufficiently low temperature 
in units of the black-hole charge, the electric field stretching to AdS-infinity causes a spontaneous 
discharge of the bulk vacuum outside of the horizon into the charged fields of the theory — whatever 
their nature. The positively charged excitations are repelled by the black hole, but cannot escape 
to infinity in AdS and they form a charge cloud hovering over the horizon. The negatively charged 
excitations fall into the black-hole and neutralize the charge, until one is left with an uncharged 
black hole with a condensate at finite T or a pure asymptotically AdS condensate solution at T = 0. 
As [2S1 l26j and we show, the statistics of the charged particle do not matter for this condensate 
formation, except in the way it forms: bosons superradiate and fermions nucleate. The dual CFT 
perspective of this process is "entropy collapse". The final state therefore has negligible ground 
state entropy and is stable. The study of charged black holes in AdS/CFT is therefore a novel way 
to understand the stability of charged interacting matter which holds much promise. 



2 From Green's function to AdS/CFT rules for a Fermi Liquid 

We wish to show how a solution with finite fermion number — a Fermi liquid — is encoded in AdS. 
The exact connection and derivation will require a review of what we have learned of Dirac field 
dynamics in AdS/CFT through Green's functions analysis. The defining signature of a Fermi liquid 
is a quasi-particle pole in the (retarded) fermion propagator, 

Gr = — + regular (2.1) 

^ — I^R — Vf{k — kp) 

Phenomenologically, a non-zero residue at the pole, Z, also known as the pole strength, is the 
indicator of a Fermi liquid state. Migdal famously related the pole strength to the occupation 
number discontinuity at the pole (w = 0): 



where 



\mi[np{uj — — npijjJ + ^)] (2.2) 



where fpo is the Fermi-Dirac distribution function. Vice versa, a Fermi liquid with a Fermi-Dirac 
jump in occupation number at the Fermi energy ujp = has a low-lying quasiparticle excitation. 
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Using our knowledge of fermionic spectral functions in AdS/CFT we shall first relate the pole- 
strength Z to known AdS quantities. Then, using Migdal's relation, the dual of a Fermi liquid is 
characterized by an asymptotically AdS solution with non-zero value for these very objects. 

The Green's functions derived in AdS / CFT are those of charged fermionic operators with scaling 
dimension A, dual to an AdS Dirac field with mass m = A — |. We shall focus on = 2 + 1 
dimensional CFTs. In its gravitational description, this Dirac field is minimally coupled to 3 + 1 
dimensional gravity and electromagnetism with the action 



S = d'^Xy/^ 



(2.3) 



For zero background fermions, \E' = 0, a spherically symmetric solution is a charged AdS4 black- hole 
background 

ds^ = —{-f{z)dt' + dx' + dy') + -j^^, 
f{z) = (l-z){l + z + z'-q'z') , 

4^3) ^ 2qaiz - 1) . (2.4) 

Here A^^^^ is the time- component of the ?7(l)-vector-potential, L is the AdS radius and the temper- 
ature and chemical potential of the black hole equal 

T=^(3-g2) , f,, = -2qa, (2.5) 

where q is the black hole charge. 

To compute the Green's functions we need to solve the Dirac equation in the background of this 
charged black hole: 



e^r^(DM + tegAM)^ + = , (2.6) 



where the vielbein , covariant derivative Dm and connection Am correspond to the fixed charged 
AdS black- hole metric and electrostatic potential (2.4) and g is the fermion charge. Denoting Aq = ^ 
and taking the standard AdS-fermion projection onto \E'-|- = ^(l±r^)\E', the Dirac equation reduces 
to 



with 



(S, + .4±)*± = tTSt (2.7) 



22 V 2/; z^' 

af a^j 

Here 7^ are the 2+1-dimensional Dirac matrices, obtained after decomposing the 3-1-1 dimensional 
F'^-matrices. 

Explicitly the Green's function is extracted from the behavior of the solution to the Dirac 
equation at the AdS-boundary. The boundary behavior of the bulk fermions is 

^^[u.k-z) = A+zt-'" + 5+^t+™ + . . . , 

^.itu,k;z) = A.zi-'" + 5_^5+™ + . . . , (2.9) 
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where A±{uj, k), B±{uj, k) are not all independent but related by the Dirac equation at the boundary 



(2m - 1) ^ + ' + (2m + 1) ^ ^ ^ 

The CFT Green's function then equals [121 ESI [H] 

= lim z-^"'^^ - singular = ^ . (2.11) 

In other words, -B_ is the CFT response to the (infinitesimal) source ^4+. Since in the Green's 
function the fermion is a fiuctuation, the functions '^±{z) are now probe solutions to the Dirac 
equation in a fixed gravitational and electrostatic background (for ease of presentation we are 
considering '^±{z) as numbers instead of two-component vectors). The boundary conditions at 
the horizon/AdS interior determine which Green's function one considers, e.g. infalling horizon 
boundary conditions yield the retarded Green's function. For non-zero chemical potential this 
fermionic Green's function can have a pole signalling the presence of a Fermi surface. This pole 
occurs precisely for a (quasi-) normalizable mode, i.e. a specific energy up and momentum kp where 
the external source A^{uj^ k) vanishes (for infalling boundary conditions at the horizon). 

Knowing that the energy of the quasinormal mode is always ujp = [TT] and following [13], we 
expand Gr around w = as: 

5(0) + ^5(1) + 

- xf^ZJ^: <^'^^' 

A crucial point is that in this expansion we are assuming that the pole will correspond to a stable 
quasiparticle, i.e. there are no fractional powers of u less than unity in the expansion around 
cui? = [13]. Fermions in AdS/CFT are of course famous for allowing more general pole-structures 
corresponding to Fermi-surfaces without stable quasiparticles [I3], but those Green's functions are 



not of the type (2.1) and we shall therefore not consider them here. The specific Fermi momentum 
kp associated with the Fermi surface is the momentum value for which the first w-independent term 
in the denominator vanishes A^^\kF) = — for this value oi k = kp the presence of a pole in the 
Green's functions at a; = is manifest. Writing A^^ = a+{k — kp) + . . . and comparing with the 
standard quasi-particle propagator, 

Gr = T- p— + regular (2.13) 

uj - fiR- vp[k - kp) 

we read off that the pole-strength equals 

Z = B^^\kp)/Al\kp). 

We thus see that a non-zero pole-strength is ensured by a non-zero value of B-{u = 0,k = kp) 
— the "response" without corresponding source as A^^\kp) = 0. Quantatively the pole-strength 
also depends on the value of A^:^\kp) = du]A^{kp)\i^=o, which is always finite. This is not a truly 
independent parameter, however. The size of the pole-strength has only a relative meaning w.r.t. to 
the integrated spectral density. This normalization of the pole strength is a global parameter rather 
than an AdS boundary issue. We now show this by proving that A^\kp) is inversely proportional 
to B^^\kp) and hence Z is completely set by B^^\kp), i.e. Z ~ \B^^\kp)\'^. Consider a transform 
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W{'^+^a,^+,b) of the Wronskian W{'^+^a,'^+,b) = '^+,Adz'^+,B - for two solutions 

to the second order equivalent of the Dirac equation for the field \1'+ 



(92 + p(^)9, + g+(^))v]/+ = o 



(2.14) 



that is conserved (detailed expressions for P{z) and Q+{z) are given in eq. (2.21)): 



W^(^+,a(^), ^+,b(^), z; zo) = exp / P(^) H^(^+,a(^), ^+,b(^)) , d,W = 0. (2.15) 



Here Zq^ is the infinitesimal distance away from the boundary at 2; = which is equivalent to 
the f/y-cutoff in the CFT. Setting k = kp and choosing for \E'+,yi = ^+2^^^""^ X]^o ^^"^ 
'^+,B = X^^o^"-^"^ ^^^^ solutions which asymptote to solutions with -B+(ci;, fci?) = 

and Aj^^oj^kp) = respectively, but for a value of u) infinitesimally away from up = 0, we can 
evaluate W at the boundary to find^ 

W = z^{l + 2m)A+B+ = iizlA+B_ (2.16) 



The last step follows from the constraint (2.10) where the reduction from two-component spinors 
to functions means that 7° is replaced by one of its eigenvalues ±i. Taking the derivative of W at 
w = for k = kp and expanding A+{u, kp) and B^{uj, kp) as in (2.12), we can solve for A^^\kp) in 



terms of B^^\kp) and arrive at the expression for the pole strength Z in terms of \B^^'{kp 



(0), 



B^^\kp)\' . 



(2.17) 



Because d^jW, as W, is a number that is independent of z, this expression emphasizes that it is truly 
the nonvanishing subleading term B^^\ujp, kp) which sets the pole strength, up to a normalization 
d^W which is set by the fully integrated spectral density. This integration is always UV-cut- 
off dependent and the explicit zq dependence should therefore not surprise us|^ We should note 
that, unlike perturbative Fermi liquid theory, Z is a dimensionful quantity of mass dimension 
2m + 1 = 2A — 2, which illustrates more directly its scaling dependence on the UV-energy scale 2:0- 
At the same time Z is real, as it can be shown that both dujW\uj=o,k=kp 
real [13]. 



/i^o'^?^-"^ and are 



4p(2) = -3/z + . . . near z = 

^Using that W is conserved, one can e.g. compute it at the horizon. There each sohition \E'+.yi(ci;, fcj?; z), 
\1/-i-_b(w, kp', z) is a linear combination of the infaUing and outgoing solution 

*+,^(z) = a(l-z)-i/^+^"/*^^ + a(l-z)-i/^-*"/*^^ + ... 
vI/+,b(z) = ^(l-z)-^/^+^"/^'^^ + /3(l-z)-i/^-^'^/^"^ + ... 

yielding a value of d^jW equal to {P{z) = 1/2(1 — z) + . . . near z = 1) 



(2.18) 



2ttT 



(2.19) 



with Af{zo) = exp dz P{z) 



2(1-2) 
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2.1 The AdS dual of a stable Fermi Liquid: Applying Migdal's relation holographi- 
cally 

We have thus seen that a solution with nonzero B_{ujp,kp) whose corresponding external source 
vanishes (by definition oi ujp, kp), is related to the presence of a quasiparticle pole in the CFT. 
Through Migdal's theorem its pole strength is related to the presence of a discontinuity of the 
occupation number, and this discontinuity is normally taken as the characteristic signature of the 
presence of a Fermi Liquid. Qualitatively we can already infer that an AdS gravity solution with 
non-vanishing B-{uf, kp) corresponds to a Fermi Liquid in the CFT. We thus seek solutions to 
the Dirac equation with vanishing external source but non-vanishing response B_ coupled to 
electromagnetism (and gravity). The construction of the AdS black hole solution with a finite single 
fermion wavefunction is thus analogous to the construction of a holographic superconductor [18] 
with the role of the scalar field now taken by a Dirac field of mass m. 

This route is complicated, however, by the spinor representation of the Dirac fields, and the 
related fermion doubling in AdS. Moreover, relativistically the fermion Green's function is a matrix 
and the pole strength Z appears in the time-component of the vector projection TiiYG. As we take 
this and the equivalent jump in occupation number to be the signifying characteristic of a Fermi 
liquid state in the CFT, it would be much more direct if we can derive an AdS radial evolution 
equation for the vector-projected Green's function and hence the occupation number discontinuity 
directly. From the AdS perspective is also more convenient to work with bilinears such as Green's 
functions, since the Dirac fields always couple pairwise to bosonic fields. 

To do so, we start again with the two decoupled second order equations equivalent to the Dirac 



equation (2.7) 



with 



{d^ + P{z)d, + Q±(z)) ^± = (2.20) 

r 



Q^iz) = + id,A±) - [d,,r\^A± + . (2.21) 

Note that both P{z) and Q±{z) are matrices in spinor space. The general solution to this second 
order equation — with the behavior at the horizon/interior appropriate for the Green's function one 
desires — is a matrix- valued function {M±{z))"ij and the field '^±{z) equals '^±{z) = M±{z)'^^^''^\ 
Due to the first order nature of the Dirac equation the horizon values ^f^j!"^^ are not independent but 
related by a z- independent matrix S'\l/^°'^'' = '$^^°^\ which can be deduced from the near-horizon 



behavior of (2.10); specifically 5* = 7 . One then obtains the Green's function from the on-shell 



boundary action (see e.g. [5^ \n\) 

Sbnd = (p (2.22) 



2 = 20 



as follows: Given a boundary source (+ for \I/+(-2), i.e. "^^Izq) = C+; one concludes that vE'^"^-' = 
M-^{zo)C+ and thus ^+{z) = M+(z)M+^(zo)C+, ^-(-2) = M_{z)SM-\zo)C+- Substituting these 
solutions into the action gives 

Sknd=(f d''xC+M4zo)SM-\zo)C+ (2.23) 

^2 = 20 
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The Green's function is obtained by differentiating w.r.t. C+ and C+ and discarding the conformal 
factor Zq"^ with m being the AdS mass of the Dirac field (one has to be careful for mL > 1/2 with 
analytic terms [34j ) 

G = lim Zo^'^M^{zo)SM-\zo) . (2.24) 

Since M±{z) are determined by evolution equations in z, it is clear that the Green's function 
itself is also determined by an evolution equation in z, i.e. there is some function G{z) which reduces 
in the limit z — )■ to Zq"^G. One obvious candidate is the function 

^(oH(^) = M_{z)SM-\z) . (2.25) 

Using the original Dirac equations one can see that this function obeys the non-linear evolution 
equation 

d,G^°^''\z) = -A^G^°^\z) - JM+SM-^ + A+G^°^''\z) + G^°^''\z)JG^°^''\z) . (2.26) 

This is the approach used in [IT], where a specific choice of momenta is chosen such that M+ 
commutes with S. For a generic choice of momenta, consistency requires that one also considers 
the evolution equation for M^{z)SM^^{z). 

There is, however, another candidate for the extension G{z) which is based on the underlying 
boundary action. Rather than extending the kernel M^{zq)M^^{zq) of the boundary action we 
extend the constituents of the action itself, based on the individual fermion wavefunctions '^±{z) = 
M±(z)S2^^ M^'^(zo). We define an extension of the matrix G{z) including an expansion in the 
complete set = {1,7*, 7*-', . . . ,7*!'*^} (with 7^ = ^7°) 

G\z) = M-\zo)M+{z)r^ M^{z)SM-\zo) , G\zo) = r^G{zo) (2.27) 

where M = i'f'^MH'j^. Using again the original Dirac equations, this function obeys the evolution 
equation 

d,G\z) = -{A+ + A_)G\z) - M-}^M_{z)fT^ M_{z)SM-}^ + M-}^M+{z)T^fM+{z)SM-}^ 

(2.28) 

Recall that fY'^-'-'-p = 7~[*i^---*p] -|- 'j~jjn---ip_ jg tj^en straightforward to see that for consistency, 
we also need to consider the evolution equations of 

Jl = M-}^M+{z)r^ M+{z)SM-}^ , Jl = M-lM_{z)r^M_{z)SM-l 

and 

G' = M-lM4z)T^ M+{z)SM-}^. 

They are 

d,ji'-''{z) = -2Re{A+)ji'-'' - r^''G''-'^\z) - TjO^'^-'^iz) - g^'^-'^-^{z)V^^ - G''-'^^ {z)rj 
d^ji'-'^iz) = -2Re(^_)X'-"^ + r^''G''-'^\z) + %G''^-'^{z) + G^'^-'^-' {z)r'^ + G'^-'^' {z)Tj 

d,&' -'-{z) = -{A- + A+)G''-'^' - r^'' j'^-''\z) - %j'^'-'\z) - j'^y-'^-\z)v-'^ + j'2-''\z)ri 

(2.29) 
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The significant advantage of these functions G^, G^, J'j. is that the evolution equations are now hn- 
ear. This approach may seem overly complicated. However, if the vector T* happens to only have a 
single component nonzero, then the system reduces drastically to the four fields J'^, G^, G^.We shall 
see below that a similar drastic reduction occurs, when we consider only spatially and temporally 
averaged functions = f dtdPxjj.. 

Importantly, the two extra currents jj. have a clear meaning in the CFT. The current G\z) 
reduces by construction to times the Green's function G^(zo) on the boundary, and clearly G\z) 
is its hermitian conjugate. The current Jl reduces at the boundary to Jl = M+^SM^]^. Thus 



J'l sets the normalization of the linear system (2.29). The interesting current is the current J7i 



Using that S = S ^, it can be seen to reduce on the boundary to the combination J'^G^T^G^. 

Thus, {J^y J} is the norm squared of the Green's function, i.e. the probability density of the 
off-shell process. 

For an off-shell process or a correlation function the norm-squared has no real functional mean- 
ing. However, we are specifically interested in solutions in the absence of an external source, i.e. 
the on-shell correlation functions. In that case the analysis is quite different. The on-shell condition 
is equivalent to choosing momenta to saturate the pole in the Green's function, i.e. it is precisely 
choosing dual AdS solutions whose leading external source A± vanishes. Then M_|_ and M_ are 
no longer independent, but M+ q = = —j^^^M_flS. As a consequence all boundary 

values of jL{zo), G^{zo), G\zo) become proportional; specifically using S" = 7° one has that 

J'^{zo)\ on- shell = ' ^'1^ G''^ {zo)\on-shell (2.30) 

/i 

is the "on-shell" Green's function. Now, the meaning of the on-shell correlation function is most 
evident in thermal backgrounds. It equals the density of states p{uj{k)) = —-^ImGji times the 
Fermi-Dirac distribution [35] 

Tm^'GUu^,Uk),k)l^_,^^,, = 2nf^n ( "^""^^^ ~ ^ ) p(c^6a.e(fc)) (2.31) 



For a Fermi liquid with the defining off-shell Green's function (2.1), we have oobareikp) — /i = w = 



and p{ujbare{k)) = Zz,j6'^{k — kF)6{oo) + . . .. Thus we see that the boundary value of J'^\zQ)\on-sheii = 
^/fd (0)5^(0) indeed captures the pole strength, times a product of distributions. This product of 
distributions can be absorbed in setting the normalization. An indication that this is correct is 
that the determining equations for G' , G^ , Jj. remain unchanged if we multiply G^, G^, Jj. on 
both sides with M+^o- If ^+,o is unitary it is just a similarity transformation. However, from the 
definition of the Green's function, one can see that this transformation precisely removes the pole. 
This ensures that we obtain finite values for G^ ^ G^, Jj. at the specific pole- values upykp where 
the distributions would naively blow up. 

2.1.1 Boundary conditions and normalizability 



We have shown that a normalizable solution to jzl^^ from the equations (2.29) correctly captures 
the pole strength directly. However, 'normalizable' is still defined in terms of an absence of a source 
for the fundamental Dirac field \E'± rather than the composite fields J'^ and G^. One would prefer 
to determine normalizability directly from the boundary behavior of the composite fields. This can 
be done. Under the assumption that the electrostatic potential $ is regular, i.e. 

<^ = fi- pz + ... (2.32) 
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the composite fields behave near z = as 



T'O _ „■ ^3~2m I „-+^4 , ■ 5+2m , 

T'O _ „■ ^5-2m I „--^4 , ■ ^3+2m , 

= Ii-2mZ'^ ^"^ + I'iZ^ + Ii+2mZ^^'^^ + -^5^^ + • ■ • , (2.33) 



with the identification 



i3-2m = |A+r, J4^ = AIB+ + 5^^+, j5+2m = \B, 
j3+2m = l^-l', JT = ^-5- + 5-^-, J5-2m = |S_ 
/4-2m = A+A^ + v4_A+, /3 = A+B^ + /4+2m = B+B^ + 

Is = B+A^ + . (2.34) 

A 'normahzable' solution in is thus defined by the vanishing of both the leading and the 
subleading term. 

3 An AdS Black hole with Dirac Hair 

Having determined a set of AdS evolution equations and boundary conditions that compute the 
pole strength Z directly through the currents J'^\z) and G\z), we can now try to construct 
the AdS dual of a system with finite fermion density, including backreaction. As we remarked 



|2 



in the beginning of section 2A, the demand that the solutions be normahzable means that the 
construction of the AdS black hole solution with a finite single fermion wavefunction is analogous 
to the construction of a holographic superconductor [18] with the role of the scalar field now taken by 



the Dirac field. The starting point therefore is the charged AdS4 black-hole background (2.4) and we 
should show that at low temperatures this AdS Reissner-Nordstrom black hole is unstable towards 
a solution with a finite Dirac profile. We shall do so in a simplified "large charge" limit where we 
ignore the gravitational dynamics, but as is well known from holographic superconductor studies 
(see e.g. [T9| [2Dj ) this limit already captures much of the essential physics. In a companion 
article [32] we will construct the full backreacted groundstate including the gravitational dynamics. 



In this large charge non-gravitational limit the equations of motion for the action (2.3) reduce 
to those of [/(l)-electrodynamics coupled to a fermion with charge g in the background of this black 
hole: 

= e^iT^{DM + iegAM)^ + mm . (3.1) 
bei n and and covariant derivative Dm remain those of the fixed charged AdS black 



Thus the vie 



hole metric (2.4), but the vector-potential now contains a background piece A^q^"* plus a first-order 



piece Am = + , which captures the effect of the charge carried by the fermions. 

Following our argument set out in previous section that it is more convenient to work with the 
currents J'].{z),G\z) instead of trying to solve the Dirac equation directly, we shall first rewrite 
this coupled non-trivial set of equations of motion in terms of the currents while at the same time 
using symmetries to reduce the complexity. Although a system at finite fermion density need not 
be homogeneous, the Fermi liquid ground state is. It therefore natural to make the ansatz that 
the final AdS solution is static and preserves translation and rotation along the boundary. As the 
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Dirac field transforms non-trivially under rotations and boosts, we cannot make tliis ansatz in tlie 
strictest sense. However, in some average sense wliicli we will make precise, the solution should be 
static and translationally invariant. Then translational and rotational invariance allow us to set 
Ai = 0, Az = 0, whose equations of motions will turn into contraints for the remaining degrees of 
freedom. Again denoting Aq = ^, the equations reduce to the following after the projection onto 



2^ 



with 



{d, + A±)^± = tT^^ (3.2) 



2zy 2fJ z^' 
r - ^^^^7° + ^^.f . (3.3) 

as before. 

The difficult part is to "impose" staticity and rotational invariance for the non-invariant spinor. 
This can be done by rephrasing the dynamics in terms of fermion current bilinears, rather than the 
fermions themselves. We shall first do so rather heuristically, and then show that the equations ob- 
tained this way are in fact the fiow equations for the Green's functions and composites J'\z), G\z) 
constructed in the previous section. In terms of the local vector currents 



J^(x, z) = ^+(a;, z)z7^*+(a;, z) , J'^lx, z) = ^_(x, z)iY'^-{x, z) , (3.4) 

or equivalently 

J^l{p^z) = j d^k^+{-k,z)ij^''^+{p + k,z) , J'^{p,z) = J d^k^_{-k, z)ij'''<i/_{p + k, z) . (3.5) 

rotational invariance means that spatial components Jj_ should vanish on the solution — this 
solves the constraint from the Ai equation of motion, and the equations can be rewritten in terms 
of only. Staticity and rotational invariance in addition demand that the bilinear momentum 
Pfj^ vanish. In other words, we are only considering temporally and spatially averaged densities: 
J±{z) = J dtd'^x^{t,x, z)i'~f^'${t,x, z). When acting with the Dirac operator on the currents to 
obtain an effective equation of motion, this averaging over the relative frequencies u and momenta 
ki also sets all terms with explicit fcj-dependence to zeroQ Restricting to such averaged currents 

^In our conventions ^ = ^£'^17". 
^ To see this consider 

{d + 2A±)'fii~k)^±{k) = Tj (*li7"*+ + *ii7"*-) + ^ - *t7'*-) ■ (3.6) 

The term proportional to $ is relevant for the solution. The dynamics of the term proportional to ki is 

{d + A++ A-){'flf'f+ - *^7'*_) = -2i^{^l-/°^+ + ¥_-f"-9^) . (3.7) 

V J 

The integral of the RHS over fc* vanishes by the assumption of translational and rotational invariance. Therefore the 



LHS of (3.7 1 and thus the second term in eq. (3.6) does so as well 
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and absorbing a factor of g/a in ^ and a factor of gy in we obtain effective equations of 
motion for the bilinears directly 



2$ 

y 



(3.8) 



with / = \E'_^+ + \E'+\&_, and all fields are real. The remaining constraint from the equation of 
motion decouples. It demands lm(iI>+^_) = - = 0. What the equations ( |3^ tell 

bulk. 



US is that for nonzero there is a charged electrostatic source for the vector potential $ in the 



Before we motivate the effective equations (3.8) at a more fundamental level, we note that these 
equations contain more information than just current conservation 9^ = 0. In an isotropic and 
static background, current conservation is trivially satisfied since d^J^ = doJ^ = —i J duje~'^^*ujJ'^{uj) 
as J^iuj 7^ 0) = 0. Also, note that we have scaled out the electromagnetic coupling. AdS4/CFT3 
duals for which the underlying string theory is known generically have g = n/L with k the grav- 



itational coupling constant as defined in (2.3). Thus, using standard AdS4/CFT3 scaling, a finite 
charge in the new units translates to a macroscopic original charge of order L/k oc A^^/^. This large 
charge demands that backreaction of the fermions in terms of its bilinear is taken into account as 
a source for 



The justification of using (3.8) to construct the AdS dual of a regular Fermi liquid is the connec- 



tion between local fermion bilinears and the CFT Green's function. The complicated flow equations 



(2.29) reduce precisely to the first two equations in (3.8) upon performing the spacetime averaging 
and the trace, i.e. J± = j d^kTiJ^ and I = J d^kTr (G^ + (5^). Combined with the demand that 



we only consider normalizable solutions and the proof that is proportional to the pole-strength. 



the radial evolution equations (3.8) are the (complicated) AdS recasting of the RG-flow for the 
pole-strength. This novel interpretation ought to dispel some of the a priori worries about our un- 
conventional treatment of the fermions through their semi-classical bilinears. There is also support 
from the gravity side, however. Recall that for conventional many-body systems and fermions in 
particular one first populates a certain set of states and then tries to compute the macroscopic 



properties of the collective. In a certain sense the equations (3.8) formulate the same program but 



in opposite order: one computes the generic wavefunction charge density with and by imposing the 
right boundary conditions, i.e normalizability, one selects only the correct set of states. This follows 
directly from the equivalence between normalizable AdS modes and quasiparticle poles that are 
characterized by well defined distinct momenta kp (for u = up = 0). The demand that any non- 
trivial Dirac hair black hole is constructed from normalizable solutions of the composite operators 
(i.e. their leading and subleading asymptotes vanishr| thus means that one imposes a superselection 



One can verify that the discussion in section 2.1.1 holds also for fully backreacted solutions. The derivation there 
builds on the assumption that the boundary behavior of the electrostatic potential is regular. It is straightforward 
to check in (3.8) that indeed precisely for normalizable solutions, i.e. in the absence of explicit fermion-sources, 
when both the leading and subleading terms in and / vanish, the boundary behavior the scalar potential remains 
regular, as required. 
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rule on the spatial averaging in the definition of Jl: 



( ^ ) I normalizable 



(fk^l ±{~k)i'^^^! ±{k)\normalizable 



(3.9) 



We see that the constraint of normalizability from the bulk point of the view, implies that one 
selects precisely the on-shell bulk fermion modes as the building blocks of the density J^. 



In turn this means that the true system that eqs. (3.8) describe is somewhat obscured by 



the spatial averaging. Clearly even a single fermion wavefunction is in truth the full set of two- 
dimensional wavefunctions whose momentum /c* has length kp- However, the averaging could just 
as well be counting more, as long as there is another set of normalizable states once the isotropic 
momentum surface \k\ = is filled. Pushing this thought to the extreme, one could even speculate 



that the system (3.8) gives the correct quantum-mechanical description of the many-body Fermi 



system: the system which gravitational reasoning suggests is the true groundstate of the charged 
AdS black hole in the presence of fermions. 

To remind us of the ambiguity introduced by spatial averaging, we shall give the boundary 
coefficient of normalizable solution for J° = / d^kj'^ a separate name. The quantity J^^{zq) is pro- 
portional to the pole strength, which via Migdal's relation quantifies the characteristic occupation 



number discontinuity atup = 0. We shall therefore call the coefficient J d^k\B_ 



\normalizable 



Anp. 



3. 0. 2 Thermodynamics 

At a very qualitative level the identification J'^\norm{z) = Anpz^'^'^"^ + . . . can be argued to follow 
from thermodynamics as well. From the free energy for an AdS dual solution to a Fermi liquid, one 
finds that the charge density directly due to the fermions is 



ptotal 



2m + 1 z, 







(3.10) 



with Zq^ the UV-cutoff as before. The cut-off dependence is a consequence of the fact that the 
system is interacting, and one cannot truly separate out the fermions as free particles. Were one to 
substitute the naive free fermion scaling dimension A = m + 3/2 = l, the cutoff dependence would 
vanish and the identification would be exact. 

We can thus state that in the interacting system there is a contribution to the charge density 
from a finite number of fermions proportional to 



PF 



—3 Anp 



+ 



(3.11) 



although this contribution formally vanishes in the limit where we send the UV-cutoff Zq ^ to infinity. 



To derive eq. (3.10), recall that the free energy is equal to minus the on-shell action of the AdS 



dual theory. Since we disregard the gravitational backreaction, the Einstein term in the AdS theory 
will not contain any relevant information and we consider the Maxwell and Dirac term only. We 
write the action as 



S 



-9 



AnDmF^^^ - ^Hp^i - m^m 



+ 



-h 



z=zo 



(3.12) 
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where we have included an exphcit fermionic boundary term that follows from the AdS/CFT dic- 
tionary [12] and is a normal vector to the boundary. The boundary action is not manifestly 
real, but its on-shell value which contributes to the free energy is real. Recall that the imaginary 
part of ^'+\E'„ decouples from eqs. (3.8). The boundary Dirac term in (3.12) is therefore equal to 
/ = 2Re('l'+^_). 

To write the free energy in terms of the quantities n, p and Anp, note that the on-shell bulk 
Dirac action vanishes. Importantly the bulk Maxwell action does contribute to the free energy. Its 
contribution is 



bulk 



lim 

2o-s>0 



/ dzd?x 




1 Zf, 


2 



- on—shell 



lim / dzd^x 

zo-s>0 



(3.13) 



on— shell 



where we have used the equation of motion (3.8). This contribution should be expected, since the 



free energy should be dominated by infrared, i.e. near horizon physics. Due to the logarithmic 



singularity in the electrostatic potential (Eq. (3.17) this bulk contribution diverges, but this diver- 
gence should be compensated by gravitational backreaction. At the same time the singularity is so 
mild, however, that the free energy, the integral of the Maxwell term, remains finite in the absence 
of the Einstein contribution. 



Formally, i.e. in the limit zq — )■ 0, the full free energy arises from this bulk contribution (3.13) 



The relation (3.10) between the charge density and Ani? follows only from the regularized free 
energy, and is therefore only a qualitative guideline. Empirically, as we will show, it is however, a 
very good one (see Fig [l] in the next section). Splitting the regularized bulk integral in two 



bulk 



dzd X 



2^377 



$(J° + J°) 



on— shell 



-\r lim 

zo->0 



dzd X 



v7 



$(4 + J°) 



(3.14) 



on— shell 



we substitute the normalizable boundary behavior of = B^z^^"^^ 
and $ = /i — + . . ., and obtain for the regularized free energy 



P Phorizoni^^*) ~l~ lim 



d xdz 



20 



2^3 



p\B_Yz' 



2^3+2m 



+ 



d^ 



X 



-B^B.zt^^^ + \ppzl 



(3.15) 



Using that 5+ = —ifi'-f^B^/{2'm + 1) (eq. (2.10)), the second bulk term and boundary contribution 
are proportional, and the free energy schematically equals 



phorizon _|_ 

20-s>0 



d^x 



3/i 



2(2m + 1 



-B_i-f°B^ 



^l+2m 
^0 



(3.16) 



With the UV-regulator Zq ^ finite, this yields the charge density in Eq. (3.10) after one recalls that 
B_ = 51^7°. 

With the derived rule that the AdS dual to a Fermi liquid has a nonzero normalizable component 
in the current J°, we will now construct an AdS solution that has this property: an AdS black hole 
with Dirac hair. Ignoring backreaction, these are solutions to the density equations (3.8). In its 
simplest form the interpretation is that of the backreaction due to a single fermion wavefunction, but 
as explained the spatial averaging of the density combined with the selection rule of normalizability 
could be capturing a more general solution. 
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3.1 At the horizon: Entropy collapse to a Lifshitz solution 



Before we can proceed with the construction of non-trivial Dirac hair solutions to Eqs. (3.8), we 



must consider the boundary conditions at the horizon necessary to solve the system. Insisting that 



the right-hand-side of the dynamical equations (3.8) is subleading at the horizon, the near-horizon 
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behavior of J^, J, $ is: 

I = 4„,(i-^)-i/2 + ...^ 

* = <:J.(1-^)M1-^) + (<1-<oI)(1-^) + ---- (3.17) 

If we insist that $ is regular at the horizon z = 1, i.e. $^^2- — O5 so that the electric field is 
finite, the leading term in must vanish as well, i.e. Jhor,± = 0, and the system reduces to a 
free Maxwell field in the presence of an AdS black hole and there is no fermion density profile 
in the bulk. Thus in order to achieve a nonzero fermion profile in the bulk, we must have an 
explicit source for the electric-field on the horizon. Strictly speaking, this invalidates our neglect of 
backreaction as the electric field and its energy density at the location of the source will be infinite. 
As we argued above, this backreaction is in fact expected to resolve the finite ground-state entropy 
problem associated with the presence of a horizon. The backreaction should remove the horizon 
completely, and the background should resemble the horizonless metrics found in [251 [371 126] ; the 
same horizon logarithmic behavior in the electrostatic potential was noted there. Nevertheless, as 
the divergence in the electric field only increases logarithmically as we approach the horizon, and 
our results shall hinge on the properties of the equations at the opposite end, near the boundary, 
we shall continue to ignore it here. We shall take the sensibility of our result after the fact as proof 
that the logarithmic divergence at the horizon is indeed mild enough to be ignored. 

The identification of the boundary value of J° with the Fermi liquid characteristic occupation 
number jump Anp rested on the insistence that the currents are built out of AdS Dirac fields. 
This deconstruction also determines a relation between the horizon boundary conditions of the 
composite fields J^, I. If '^±{z) = C±{1 - z)'^!'^ + . . . then Jhor,± = C\ and I^or = C+C^. As 
the solution is independent of the solution $1^^^^ which is regular at the horizon, we match the 
latter to the vector-potential of the charged AdS black hole: = ~2gq = gfio/a. Recalling 

that = ~{Jhor,+ + Jhor,~), we See that the three-parameter family of solutions at the horizon 

(2) 

in terms of C±, corresponds to the three-parameter space of boundary values i?_ and /i 
encoding a fermion-source, the fermion-response/expectation value and the chemical potential. 

We can now search whether within this three-parameter family a finite normalizable fermion 
density solution with vanishing source = exists for a given temperature T of the black hole. 



3.2 A BH with Dirac hair 

The equations are readily solved numerically with a shooting method from the horizon. We consider 
both an uncharged AdS-Schwarzschild solution and the charged AdS Reissner-Nordstrom solution. 
Studies of bosonic condensates in AdS/CFT without backreaction have mostly been done in the 
AdS-Schwarzschild (AdSS) background ([IHlIin] and references therein). An exception is [38j, which 
also considers the charged RN black hole. As is explained in [3S], they correspond to two different 
limits of the exact solution: the AdSS case requires that Anp ^ /i that is, the total charge of the 
matter fields should be dominant compared to the charge of the black hole. On the other hand, the 
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RN limit is appropriate if ^np <^ /x. It ignores the effect of the energy density of the charged matter 
sector on the charged black hole geometry. The AdS Schwarzschild background is only reliable near 
Tc, as at low temperatures the finite charged fermion density is comparable to /x. The RN case is 
under better control for low temperatures, because near T = the chemical potential can be tuned 
to stay larger than fermion density. 

We shall therefore focus primarily on the solution in the background of an AdS RN black hole, 
i.e. the system with a heat bath with chemical potential a — non-linearly determined by the value 

('2') 

of = /io at the horizon — which for low T//i should show the characteristic ^np of a Fermi 
liquid. The limit in which we may confidently ignore backreaction is /^o for T < — for 

AdSS the appropriate limit is ^^q*^ <^ T for /io ^ T. 



3.2.1 Finite fermion density solutions in AdS-RN 

Fig. [l] shows the behavior of the occupation number discontinuity up = l-B-P and the fermion 
free-energy contribution / as a function of temperature in a search for normalizable solutions to 
Eqs ( |3.8[ ) with the aforementioned boundary conditions. We clearly see a first order transition to 
a finite fermion density, as expected. The underlying Dirac field dynamics can be recognized in 
that the normalizable solution for J^{z) which has no leading component near the boundary by 
construction, also has its subleading component vanishing (Fig. [2| 
Analyzing the transition in more detail in Fig. [3} we find: 

1. The dimensionless number discontinuity Anp/fi'^'^ scales as in a certain temperature 
range Tp < T < Tc, with 6 > depending on g and A, and Tp typically very small. At 
T = Tc > Tp it drops to zero discontinuously, characteristic of a first order phase transition. 

2. At low temperatures, < T < Tp, the power-law growth comes to a halt and ends with a 
plateau where Anp/ fi^^ ~ const. (Fig. |3]A). It is natural to interpret this temperature as the 
Fermi temperature of the boundary Fermi liquid. 

3. The fermion free energy contribution I / jj?^^^ scales as T^/^ with u > 1 for < T < Tc, and 
drops to zero discontinuously at T^. As / empirically equals minus the free energy per particle, 
it is natural that I{T = 0) = 0, and this in turn supports the identification of Anp{T = 0) as 
the step in number density at the Fermi energy. 

One expects that the exponents 6, v are controlled by the conformal dimension A,^*^ The depen- 
dence of the exponent 5 on the conformal dimension is shown in Fig. |3]A.. While a correlation clearly 
exists, the data are not conclusive enough to determine the relation b = 5(A). The clean power law 
T~^ scaling regime is actually somewhat puzzling. These values of the temperature, Tp < T < T^ 



Cl 



correspond to a crossover between the true Fermi liquid regime for T < Tp and the conformal phase 
for T > Tc, hence there is no clear ground for a universal scaling relation for 6, which seems to 
be corroborated by the data (Fig. [3|3). At the same time, the scaling exponent u appears to obey 
u = 2 with great precision (Fig. |3b, inset) independent of A and g. 



^Although the Dirac hair solution has charged matter in the bulk, there is no Higgs effect for the bulk gauge 
field, and thus there is no direct spontaneous symmetry breaking in the boundary. Indeed one would not expect it 
for the Fermi liquid groundstate. There will be indirect effect on the conductivity similar to [IS]- We thank Andy 
O'Bannon for his persistent inquiries to this point. 

^'^The charge g of the underlying conformal fermionic operator scales out of the solution. 
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Figure 1: (A) Temperature dependence of the Fermi liquid occupation number discontinuity An^ and 
operator I for a fermionic field of mass m = —1/4 dual to an operator of dimension A = 5/4. We see 
a large density for r//i small and discontinuously drop to zero at T ~ 0.05//. At this same temperature, 
the proxy free energy contribution per particle (the negative of /) vanishes. (B) The free energy F = 
pfermion _j_ pMaxwell (^(^^ (3.12)) as a function of T//! ignoring the contribution from the gravitational 
sector. The blue curve shows the total free energy F = pMaxweii^ which is the sum of a bulk and a 
boundary term. The explicit fermion contribution Ffgrmion vanishes, but the effect of a non-zero fermion 
density is directly encoded in a non-zero F^I^^^^^^K The figure also shows this bulk F^l^^^^^ and the 
boundary contribution F^^l^^^^^ separately and how they sum to a continuous Ftot^i- Although formally 



the explicit fermion contribution i^j ~ / in equation (3.16) vanishes, the bulk Maxwell contribution is 



captured remarkably well by its value when the cut-off is kept finite. The light-green curve in the figure 
shows Ff for a finite zq ~ 10~^. For completeness we also show the total charge density, Eq. (3.10). 
dimension of the fermionic operator used in this figure is A = 1.1. 



The 



A final consideration, needed to verify the existence of a finite fermion density AdS solution 
dual to a Fermi liquid, is to show that the ignored backreaction stays small. In particular, the 
divergence of the electric field at the horizon should not affect the result. The total bulk electric 
field Ez = —dz^ is shown in Fig. |4jA., normalized by its value at z = 1/2. The logarithmic 
singularity at the horizon is clearly visible. At the same time, the contribution to the total electric 
field from the charged fermions is negligible even very close to the horizon p] This suggests that 
our results are robust with respect to the details of the IR divergence of the electric field. 

The diverging backreaction at the horizon is in fact the gravity interpretation of the first order 
transition at T^: an arbitrarily small non-zero density leads to an abrupt change in the on shell bulk 
action. As the latter is the free energy in the CFT, it must reflect the discontinuity of a first order 
transition. A full account of the singular behavior at the horizon requires self-consistent treatment 
including the Einstein equations. At this level, we can conclude that the divergent energy density at 
the horizon implies that the near-horizon physics becomes substantially different from the AdS2 limit 
of the RN metric. It is natural to guess that the RN horizon disappears completely, corresponding 

-'^^It is of the order lO"**, starting from z — 0.9999. We have run our numerics using values between 1 — 10^^ and 
1 — 10"^ and found no detectable difference in quantities at the boundary. 
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Figure 2: The boundary behaviour of J_ (0) in for a generic solution (blue) to Eqs. ( 3.8 ) and a normalizable 
Dirac-hair solution (red) for m = —1/4 in the background of an AdS-RN black hole with /i/T = 128.8. 
The dotted lines show the scaling z^^/^ and of the leading and subleading terms in an expansion of 
j£(z) near z = 0; the dashed line shows the scaling z^/^ of the subsubleading expansion whose coefficient 
is \B^{iOF, ^f)P- That the Dirac hair solution (red) scales as the subsubleading solution indicates that the 
current j£ faithfully captures the density of the underlying normalizable Dirac field. 



to a ground state with zero entropy, as hypothesized in [25] . This matches the expectation that the 
finite fermi-density solution in the bulk describes the Fermi-liquid. The underlying assumption in 
the above reasoning is that the total charge is conserved. 

3.2.2 Finite fermion density in AdSS 

For completeness, we will describe the finite fermion-density solutions in the AdS Schwarzschild 
geometry as well. In these solutions the charge density is set by the density of fermions alone. 
They are therefore not reliable at very low temperatures T ^ Tc when gravitational backreaction 
becomes important. The purpose of this section is to show the existence of finite density solutions 
does not depend on the presence of a charged black-hole set by the horizon value ^^^^^ = /zq, but 
that the transition to a finite fermion density can be driven by the charged fermions themselves. 

Fig. |4j3 shows the nearly instantaneous development of a non- vanishing expectation value for the 
occupation number discontinuity Anp and I in the AdS Schwarzschild background. The rise is not 
as sharp as in the RN background. It is, however, steeper than exponential, and we may conclude 
that the system undergoes a discontinuous first order transition to a AdS Dirac hair solution. The 
constant limit reached by the fermion density as T — )■ has no meaning as we cannot trust the 
solution far away from Tc. 

The backreaction due to the electric field divergence at the horizon can be neglected, for the 
same reason as before (Fig. |4p). 

3.3 Confirmation from fermion spectral functions 

If, as we surmised, the finite fermion density phase is the true Fermi-liquid-like ground state, the 
change in the fermion spectral functions should be minimal as the characteristic quasi-particle peaks 
are already present in the probe limit, i.e. pure AdS Reissner- Nordstrom pTl [12]. Fig. [s] shows 
that quasiparticle poles near u = with similar analytic properties can be identified in both the 
probe limit, i.e. pure AdS-RN case and the AdS-RN Dirac-hair solution. The explanation for this 
similarity is that the electrostatic potential $ almost completely determines the spectrum, and the 
change in $ due to the presence of a finite fermion density is quite small. Still, one expects that 
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Figure 3: (A) Approximate power-law scaling of the Fermi liquid characteristic occupation number 
discontinuity Anp/ii^^ ~ T^^ as a function of T/ for A = 5/4. This figure clearly shows the saturation 
of the density at very low T/fi. The saturation effect is naturally interpreted as the influence of the 
characteristic Fermi energy. (B) The scaling exponent 6 for different values of the conformal dimension 
A. There is a clear correlation, but the precise relation cannot be determined numerically. The scaling 
exponent of the current I / ji'^^'^^ ~ j^-i/u obeys v = 2 with great accuracy, on the other hand (Inset). 



the finite fermion density system is a more favorable state. This indeed follows from a detailed 
comparison between the spectral functions ^4(0;; k) in the probe limit and the fermion-liquid phase 
(Fig. g. We see that: 

1. All quasiparticle poles present in the probe limit are also present in the Dirac hair phase, at a 
slightly shifted value of kp. This shift is a consequence of the change in the bulk electrostatic 
potential $ due to the presence of the charged matter. For a Fermi-liquid-like quasipar- 
ticle corresponding to the second pole in the operator with A = 5/4 and g = 2 we find 

non-Fermi-liquid pole, i.e. the first pole for the same conformal 
operator, has k^°'"^ — fc^"^ = 0.03/i. 

2. The dispersion exponents u defined through [oj — EfY' ~ (k — kpY^'^, also maintain roughly the 
same values as both solutions. This is visually evident in the near similar slopes of the ridges 
in Fig. |5j In the AdS Reissner-Nordstrom background, the dispersion coefficients are known 

analytically as a function of the Fermi momentum: z/^^ = y^2^ — | + | (A — 3/2)^ [13] . The 
Fermi-liquid-like quasiparticle corresponding to the second pole in the operator with A = 5/4 
and g = 2 has I'l^J'^'^ = 1.02 vs. u^^^ = 1.01. The non-Fermi-liquid pole corresponding to the 
first pole for the same conformal operator, has z^^™^'^ ~ 0.10, and z/^"^ = 0.12. 

3. The most distinct property of the finite density phase is the redistributed spectral weight 
of the poles. The non-Fermi liquid pole reaches its maximum height about 10^, an order of 
magnitude less than in the probe limit, whereas the second, Fermi liquid-like pole, increases 
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Figure 4: (A) The radial electric field —Ez = d^/dz, normalized to the midpoint value Ez{z)/Ez{l/2) 
for whole interior of the finite fermion density AdS-RN solution (upper) and near the horizon (lower). 
One clearly sees the soft, log-singularity at the horizon. The colors correspond to increasing temperatures 
from T = 0.04/i (lighter) to T = 0.18/u (darker), all with A = 1.1. (B) The occupation number jump 
Anp and free energy contribution / as a function of temperature in AdS-Schwarzschild. We see the jump 
Anp saturate at low temperatures and fall off at high T. An exponential fit to the data (red curve) 
shows that in the critical region the fall-off is stronger than exponential, indicating that the transition 
is first order. The conformal dimension of the fermionic operator is A = 1.1. (C) The radial electric 
field —Ez = d^/dz, normalized to the midpoint value {Ez{z)/Ez{l/2)) for the finite fermion density AdS- 
Schwarzschild background. The divergence of the electric field Ez is again only noticeable near the horizon 
and can be neglected in most of the bulk region. 



by an order of magnitude. This suggests that the finite density state corresponds to the 
Fermi-hquid hke state, rather than a non- Fermi hquid. 

4. As we mentioned in the introduction, part of the reason to suspect the existence of an AdS- 
RN Dirac-hair solution is that a detailed study of spectral functions in AdS-RN reveals that 
the quasiparticle peak is anomalously sensitive to changes in T. This anomalous temperature 
dependence disappears in the finite density solution. Specifically in pure AdS-RN the position 
oJmax where the peak height is maximum, denoted Ep in [12], does not agree with the value 
(^poie, where the pole touches the real axis in the complex cj-plane, for any finite value of T, 
and is exponentially sensitive to changes in T (Fig [6]). In the AdS-RN Dirac hair solution the 
location Umax and the location Upoie do become the same. Fig. [6p shows that the maximum 
of the quasiparticle peak always sits at w ~ in finite density Dirac hair solution, while it 
only reaches this as T — )■ in the probe AdS-RN case. 

4 Discussion and Conclusion 

Empirically we know that the Fermi liquid phase of real matter systems is remarkably robust and 
generic. This is corroborated by analyzing effective field theory around the Fermi surface, but as 
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Figure 5: The single-fermion spectral function in the probe hmit of pure AdS Reissner-Nordstrom 
(red/yeUow) minus the spectrum in the finite density system (blue). The conformal dimension is A = 5/4, 
the probe charge g = 2, and fi/T = 135. We can see two quasiparticle poles near a; = 0, a non-FL pole 
with A:^™'"' ~ 0.11/i and /c^"^ ~ 0.08/i respectively and a FL-pole with k^p"^'' ~ 0.18/i and fc^"-^ ~ O.Ufj,. 
The dispersion of both poles is visibly similar between the probe and the finite density backgroudnd. At 
the same time, the non-FL pole has about 8 times less weight in the finite density background, whereas 
the FL-pole has gained about 6.5 times more weight. 



it assumes the ground state it cannot explain its genericity. If the Fermi hquid ground state is 
so robust, this must also be a feature of the recent holographic approaches to strongly interacting 
fermionic systems. Our results here indicate that this is so: We have used Migdal's relation to 
construct AdS/CFT rules for the holographic dual of a Fermi liquid: the characteristic occupation 
number discontinuity Anp is encoded in the normalizable subsubleading component of the spatially 
averaged fermion density Jl{z) = J (fk'^{u = 0, —k, z)i'j'^'^{u = 0,k,z) near the AdS boundary. 
This density has its own set of evolution equations, based on the underlying Dirac field, and insisting 
on normalizability automatically selects the on-shell wavefunctions of the underlying Dirac-field. 

The simplest AdS solution that has a non- vanishing expectation value for the occupation number 
discontinuity Anp is that of a single fermion wavefunction. Using the density approach — which 
through the averaging appears to describe a whole class of solutions rather than one specific solution 
— we have constructed the limit of this solution where gravitational backreaction is ignored. At 
low black hole temperatures this solution with fermionic "Dirac hair" is the preferred ground state. 
Through an analysis of the free-energy, we argue that this gravitational solution with a non-zero 
fermion profile precisely corresponds to a system with a finite density of fermions. A spectral anal- 
ysis still reveals a zoo of Fermi-surfaces in this ground state, but there are indications that in the 
full gravitationally backreacted solution only a Landau Fermi-liquid type Fermi surface survives. 
This follows in part from the relation between the spectral density and the Fermi momentum of a 
particular Landau liquid-like Fermi surface; it also agrees with the prediction from Luttinger's theo- 
rem. Furthermore, the spectral analysis in the finite density state shows no anomalous temperature 
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Figure 6: (A) Single fermion spectral functions near a; = in pure AdS Reissner-Nordstrom (blue) and 
in the finite fermion density background (red). In the former the position of the maximum approaches 
a; = as T is lowered whereas in the latter the position of the maximum stays close to T = for all values 
of T. (B) Position of the maximum of the quasiparticle peak in k-co plane, for different temperatures and 
A = 5/4. The probe limit around a AdS-RN black hole (blue) carries a strong temperature dependence of 
the ^max value, with oJmax,T^o 7^ 0. In the finite fermion density background, the position of the maximum 
(red) is nearly independent of temperature and stays at cj = 0. 



dependence present in the pure charged black-hole single spectral functions. This also indicates 
that the finite density state is the true ground state. 

The discovery of this state reveals a new essential component in the study of strongly coupled 
fermionic systems through gravitational duals, where one should take into account the expectation 
values of fermion bilinears. Technically the construction of the full gravitationally backreacted 
solution is a first point that is needed to complete our finding. This is under current investigation. 
The realization, however, that expectation values of fermion bilinears can be captured in holographic 
duals and naturally encode phase separations in strongly coupled fermion systems should find a large 
set of applications in the near future. 
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